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Abstract 



The quark-meson coupling model for nuclear matter, which describes nuclear 
matter as non-overlapping MIT bags bound by the self-consistent exchange 
of scalar and vector mesons, is modified by introducing medium modification 
of the bag constant. We model the density dependence of the bag constant in 
two different ways: one invokes a direct coupling of the bag constant to the 
scalar meson field, and the other relates the bag constant to the in-medium 
nucleon mass. Both models feature a decreasing bag constant with increas- 
ing density. We find that when the bag constant is significantly reduced in 
nuclear medium with respect to its free-space value, large canceling isoscalar 
Lorentz scalar and vector potentials for the nucleon in nuclear matter emerge 
naturally. Such potentials are comparable to those suggested by relativistic 
nuclear phenomenology and finite-density QCD sum rules. This suggests that 
the reduction of bag constant in nuclear medium may play an important role 
in low- and medium-energy nuclear physics. 
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I. INTRODUCTION 



Ultimately, the physics of nuclear matter and finite nuclei is an exercise in applied quan- 
tum chromodynamics (QCD), which governs the underlying strong interactions of quarks 
and gluons. In reality, however, knowledge of QCD has had very little impact, to date, 
on the study of low- and medium-energy nuclear phenomena. The reason is that QCD is 
intractable at the nuclear physics energy scales due to the nonperturbative features of QCD. 
A reasonable consensus is that the relevant degrees of freedom for low energy QCD are 
hadrons instead of quarks and gluons. 

While the description of nuclear phenomena has been efficiently formulated using the 
hadronic degrees of freedom, new challenges arise from the observed small but interest- 
ing corrections to the standard hadronic picture such as the EMC effect which reveals the 
medium modification of the internal structure of nucleon To address these new chal- 
lenges, it is necessary to build theories that incorporate quark-gluon degrees of freedom, yet 
respect the established theories based on hadronic degrees of freedom. 

A few years ago, Guichon || proposed a quark-meson coupling (QMC) model to inves- 
tigate the direct "quark effects" in nuclei. This model describes nuclear matter as non- 
overlapping MIT bags interacting through the self-consistent exchange of mesons in the 
mean-field approximation. This simple model was refined later by including nucleon Fermi 
motion and the center-of-mass corrections to the bag energy and applied to variety of 
problems |§-[7|. Recently, the QMC model has been applied to finite nuclei [§,0. (There 
have been several works that also discuss the quark effects in nuclei, based on other effective 
models for the nucleon flO| ) . 

Although it provides a simple and attractive framework to incorporate the quark struc- 
ture of the nucleon in the study of nuclear phenomena, the QMC model has a serious 
short-coming. It predicts much smaller scalar and vector potentials for the nucleon than 
obtained in relativistic nuclear phenomenology. Unless there is a large isoscalar anomalous 
coupling (ruled out by other considerations) this implies a much smaller nucleon spin-orbit 
force in finite nuclei. To lowest order in the nucleon velocity and potential depth the nu- 
cleon spin-orbit potential can be obtained in a model independent way from the strengths 
of the scalar and vector potentials. The spin-orbit potential from the QMC model is too 
weak to successfully explain spin-orbit splittings in finite nuclei and the spin-observables in 
nucleon-nucleus scattering. 

Relativistic nuclear phenomenology is a general approach based on nucleons and mesons 
and has gained tremendous credibility during last twenty years. In this framework, the nu- 
cleons in nuclear environment are treated as point-like Dirac particles interacting with large 
canceling isoscalar Lorentz scalar and vector potentials. This approach has been successful 
in describing the spin-observables of nucleon-nucleus scattering in the context of relativistic 
optical potentials JTT| , [T!2| ]. Moreover, such potentials can be derived from the relativistic 



impulse approximation ||12|| . The relativistic field-theoretical models based on nucleons and 
mesons, QHD, also feature Dirac nucleons interacting through the exchange of scalar and 
vector mesons fl3[|. QHD, at the mean-field level, has proven to be a powerful tool for de- 



scribing the bulk properties of nuclear matter and spin-orbit splittings of finite nuclei |13[ . It 



is known that the large and canceling scalar and vector potentials are central to the success 
of the relativistic nuclear phenomenology. Recent progress in understanding the origin of 
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these large potentials for propagating nucleons in nuclear matter has been made via the 



analysis of the finite-density QCD sum rules |14[ . 

In a recent paper [TjJ, the present authors have pointed out that the resulting small 
nucleon potentials in the QMC model stem from the assumption of fixing the bag constant 
at its free-space value, and that this assumption is questionable. We then included a density 
dependent bag constant and found that when the bag constant drops significantly in nuclear 
matter relative to its free-space value, the large potentials for nucleons in nuclear matter, 
as seen in the relativistic nuclear phenomenology and finite-density QCD sum rules, can be 
recovered. This suggests that the reduction of bag constant in nuclear matter relative to its 
free-space value may be essential for the successes of relativistic nuclear phenomenology and 
thus may play an important role in low- and medium-energy nuclear physics. In the present 
paper, we present further details and model the density dependence of the bag constant in 
two different ways: one invokes a direct coupling of the bag constant to the scalar meson 
field, and the other relates the bag constant to the in-medium nucleon mass. 

This paper is organized as follows: In Section [H| we sketch the QMC model for nuclear 
matter. We then modify the QMC model by introducing a density dependent bag constant 
in nuclear matter in Section |IHl The results are presented in Section [TV]. Further discussions 
are given in Section |VJ Section [VI| is a summary. 



II. THE QUARK-MESON COUPLING MODEL FOR NUCLEAR MATTER 



In this section, we give a brief introduction to the quark-meson coupling model for nuclear 
matter. The reader is referred to Refs. [@-[5| for further details and justifications for using a 
simple bag model for the in-medium nucleon. 

In the QMC model, the nucleon in nuclear medium is assumed to be a static spherical 
MIT bag in which quarks interact with the scalar and vector fields, a and uJ, and these 
fields are treated as classical fields in the mean field approximation. (Here we only consider 



up and down quarks.) 
motion: 



The quark field, i/) q (x), inside the bag then satisfies the equation of 



i ft-(m q - g q c a) 



%{x) = 



(2.1) 



where m q is the current quark mass, and g q a and g% denote the quark-meson coupling con- 
stants. We will neglect isospin breaking and take m" = (mP u + m®)/2 hereafter. The nor- 
malized ground state for a quark in the bag is given by 



jo(xr/R) 
i(5 q a ■ rj^xr/R) 



X q 



4tt 



(2.2) 



where 



e q = tt q + gluR 



(3 q 



£l q — Rm* 
\\n q + Rm* 



2R 3 j%(x) n g {Q q - 1) + Rml/2 / 



'x 



(2.3) 
(2.4) 
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with Q q = ^jx 2 + (Rm*) 2 , m* = m° — g%a, R the bag radius, and Xq the quark spinor. The 
x value is determined by the boundary condition at the bag surface 

3o{x)=MM ■ (2-5) 
The energy of a static bag consisting of three ground state quarks can be expressed as 



^ba g = 3^-| + ^7T J R 3 J B, (2.6) 



where Z is a parameter which accounts for zero-point motion and B is the bag constant. In 
the discussions to follow, we use Rq, B and Z to denote the corresponding bag parameters 
for the free nucleon. After the corrections of spurious center-of-mass motion in the bag, the 
effective mass of a nucleon bag at rest is taken to be MM 



M* N = vX g - (PL) > (2-7) 



where (pl m ) = J2 q (Pq) an d (p 2 ) is the expectation value of the quark momentum squared, 
{x/Rf. 

The equilibrium condition for the bag is obtained by minimizing the effective mass 
with respect to the bag radius 



dM* N 
dR 



. (2.8) 



In free space, one may fix at its experimental value 939 MeV and use the equilibrium con- 
dition to determine the bag parameters. For several choices of bag radius, R = 0.6, 0.8, 1.0 
fm, the results for S 1/4 and Z are 188.1, 157.5, 136.3 MeV and 2.030, 1.628, 1.153, respec- 
tively. 

The total energy per nucleon at finite density p^, including the Fermi motion of the 
nucleons, can be written as 01 



Eu* = t F d 3 k^M* N 2 + k 2 + ^p N + ^U 2 , (2.9) 

(27ryp N J v 2mi 2 p N 

where 7 is the spin-isospin degeneracy, and 7 = 4 for symmetric nuclear matter and 7 = 2 
for neutron matter.^ Note that the mean field u; created by uniformly distributed nucleons 
is determined by baryon number conservation to be 



^=^% = ^, (2.10) 
where g w = 3g^. The scalar mean field is determined by the thermodynamic condition 



1 Here we only consider symmetric nuclear matter and neutron matter. The generalization to a 
general asymmetric nuclear matter is straightforward (see, for example, Ref. B). 
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'dR 



tot 



da 







R,PN 



If one assumes 

B = Bq 

and Z = Z , Eq. fl2.11|) yields the self-consistency condition 



— 9 a ^ i — \ T 
mr a [lit) 6 



kp 



d 3 k 



M 



N 



M* N 2 + k 2 



with 

9a C (a) 
where g a = 3 g% and 



E 



9a 



bag 



M* N 



1 - 



S(a) + 



m' 



Eh 



bag 



S(w) 



n q /2 + Rm* q (Q q -l) 
Vl q ((],-l) + i?m?/2 



(2.11) 



(2.12) 



(2.13) 



(2-14) 



(2.15) 



The two coupling constants g a and g u can be chosen to fit the nuclear matter binding energy 
at the saturation density. For a given density, Eqs. fl2.5Q, ( |2.8| ), and ( |2.13| ) form a set of 
equations for calculating x, R, and a. 



III. MODIFIED QUARK-MESON COUPLING MODEL 

In this section, we modify the QMC model by introducing a density dependent bag 
constant. We propose two models for the modification of the bag constant, featuring a 
decreasing bag constant with increasing nuclear matter density. In principle, the parameter 
Z may also be modified in the nuclear medium. However, it is unclear how Z changes 
with the density. Here we assume that the medium modification of Z is small at low and 
moderate densities and take Z = Z .Q 



A. Direct coupling model 

The bag constant in the MIT bag model contributes ~ 200 — 300 MeV to the nucleon 
energy and provides the necessary pressure to confine the quarks. Thus, the bag constant 
is an inseparable ingredient of the bag picture of a nucleon. When a nucleon bag is put into 



2 Recently, Blunden and Miller || have considered a density dependent Z. However, it is found 
that for reasonable parameter ranges changing Z has little effect and tend to make the model 
worse. 
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the nuclear medium, the bag as a whole reacts to the environment. As a result, the bag 
constant may be modified. There is little doubt that at sufficiently high densities, the bag 
constant is eventually melted away and quarks and gluons become the appropriate degrees 
of freedom. Therefore, It is reasonable to believe that the bag constant is modified and 
decreases as density increases. This physics is obviously bypassed in the QMC model by the 
assumption of B = B . 

To reflect this physics, we modify the QMC model by introducing a direct coupling 
between the bag constant and the scalar mean field 



B 
B~ 



B 



1- 9t 



4 a 



5M 



N 



(3.1) 



where and 5 are real positive parameters and the introduction of Mjy is based on the 
consideration of dimension. (The case, 5 — 1, was also considered by Blunden and Miller 
||.) Note that g„ differs from g% (or g a ). When g„ = 0, the usual QMC model is recovered. 

This direct coupling can be partially motivated from considering a non-topological soli- 
ton model for the nucleon where a scalar soliton field provides the confinement of the quarks. 
Roughly speaking, the bag constant in the MIT bag model mimics the effect of the scalar 
soliton field in the soliton model. Now, when a nucleon soliton is put into nuclear envi- 
ronment, the scalar soliton field will interact with the scalar mean field (see, for example, 
Refs. [16,|10[). Therefore, it is reasonable to couple the bag constant directly to the scalar 
mean field. 

The factor C(a) of (|2.14|) , appearing in the self-consistency condition Eq. (|2.13|) , then 
becomes 



g*C{o) = g, 



bag 



M* N 



fir, 



Eba. K R . 



S(W) 



ml 



+ 9, 



B Eba,, 



M 



16 D 3 B 
71 R 6 



N 



M, 



N 



l 4 gftf 

5 M; 



N 



(3.2) 



The other equations are not affected. In the limit of 5 
exponential form with a single parameter g^ 



B 
B~ n 



-4gj? a/M N 



oo, Eq. ( |3.ip reduces to an 



(3.3) 



In the limit of zero current quark mass (i.e., m = 0) and g a = 0, the nucleon mass scales 
like S 1 / 4 from dimensional arguments (see also Appendix A). Then from Eq. ( |3.1| ) we get 

M* N /M N = (B/B ) 



1/4 



1 



y°- s m n 



5/4 



We observe that the linear a-nucleon coupling 

is just g„ while 5 controls the non-linearities. For 5 = 4, the non-linearities vanish and, as 
discussed in the next section, we recover QHD-I but with a density dependent bag radius. 



B. Scaling model 



In the previous paper [15], we have considered a scaling model, which relates the in- 



medium bag constant to the in-medium nucleon mass directly 
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B _ rM 
~B~v ~ [M 

where k is a real positive parameter and k 
factor C(a), in this case, is given by 



N 



N- 



(3.4) 



corresponds to the usual QMC model. The 



g a C(a) = g a 



bag 



M 



N 



a, 



-Ew R . 



S(a) + 



m' 



1 



E] 



K 



bag 



M$ 3 



ttR 3 B 



(3.5) 



7 bag J 

Note that in this model, the effective nucleon mass and the bag constant B are deter- 
mined self-consistently by combining Eqs. (|2.6|), (|2.7fl , and ( [3.4Q . 

One notices that both Eqs. (|3.1|) and ( ft.4| ) give rise to a reduction of the bag constant 
in nuclear medium relative to its free-space value. While the scaling model is characterized 
by a single free parameter k, it leads to a complicated and implicit relation between the 
bag constant and the scalar mean field. On the other hand, the direct coupling model 
features a straightforward coupling between the bag constant and the scalar mean field, 
which, however, introduces two free parameters, g^ and S. 



IV. RESULTS 

In this section, we present numerical results. The two models for the in-medium bag 
constant discussed in previous section will be considered. The current quark masses are 
taken to be m u = = for simplicity. Inclusion of small current quark masses only leads 
to numerically small refinement of present results. 

Let us start from the direct coupling model. For a given value of g%, we adjust the 
coupling constants g„ and g u to reproduce the nuclear matter binding energy (—16 MeV) 
at the saturation density (p^, =0.17 fm~ 3 ). The resulting coupling constants and nuclear 
matter results are given in Table | for various 5 values with Rq = 0.6 fm. For the special 
case, g q a = and 5 = 4, the present model leads to exactly the same nuclear matter results 
as obtained in QHD-I (see first row of Table |). This is also shown analytically in Appendix 
A. 

The most important feature is that the reduction of B relative to B Q leads to the decrease 
of M]^/Mjv and the increase of U v /M^ relative to their values in the simple QMC model. 
In the usual QMC model, the required vector coupling is very small. This, in Refs. is 
attributed to the repulsion provided by the center-of-mass corrections to the bag energy. In 
our modified QMC model, the reduction of the bag constant in nuclear medium provides a 
new source of attraction as it effectively reduces Mjy-. Consequently, additional vector field 
strength is required to obtain the correct saturation properties of nuclear matter. 

The above physics is clearly reflected in Table ID. It can be seen from Table | that when 
the bag constant is reduced significantly in nuclear matter relative to its free-space value, 
the resulting magnitudes for — and U v = gjd are qualitatively different from those 
obtained in the simple QMC model. In particular, for 5 = (13.0, 8.0, 3.6), corresponding to 
gl = (0, 1.0, 2.0), we get B/B ~ 0.36 and 

M* N ~ 660 - 720 MeV , (4.1) 
U v ~ 150 -215 MeV , (4.2) 
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TABLE I. Coupling constants and nuclear matter results as obtained from the direct coupling 
model. The free space bag radius is fixed at Rq = 0.6 fm. Here the nuclear matter compressibility, 
K y , is given in unit of MeV, r* n and r m denote the quark root-mean-square radius in nuclear 
matter and in free space, respectively, and U v = g^uJ is the vector mean field. The mass parameters 
are taken to be m q = 0, m a = 550 MeV, and = 783 MeV. 



9% 


5 






M* N /M N 


U V /M N 




B/B 


x/xq 


R/Rq 


r m/ r m 





4 


8.45 


12.84 


0.55 


0.37 


540 


0.09 


1.0 


1.83 


1.83 




8 


5.68 


6.46 


0.75 


0.18 


313 


0.31 


1.0 


1.34 


1.34 




12 


5.40 


5.68 


0.77 


0.16 


295 


0.35 


1.0 


1.30 


1.30 




13 


5.28 


5.57 


0.77 


0.16 


293 


0.36 


1.0 


1.29 


1.29 




oc 


4.95 


4.62 


0.80 


0.13 


270 


0.41 


1.0 


1.25 


1.25 


1.0 


4 


5.69 


10.84 


0.61 


0.32 


490 


0.19 


0.97 


1.51 


1.52 




8 


4.20 


6.78 


0.74 


0.19 


333 


0.36 


0.97 


1.28 


1.29 




12 


3.96 


6.14 


0.76 


0.18 


315 


0.39 


0.97 


1.25 


1.26 




oc 


3.69 


5.24 


0.78 


0.15 


289 


0.45 


0.98 


1.22 


1.22 


2.0 


3.6 


3.16 


8.03 


0.70 


0.23 


431 


0.36 


0.93 


1.27 


1.29 




4 


2.99 


7.42 


0.72 


0.21 


398 


0.39 


0.94 


1.25 


1.27 




8 


2.54 


5.81 


0.77 


0.17 


336 


0.48 


0.95 


1.18 


1.20 




12 


2.43 


5.48 


0.76 


0.16 


324 


0.50 


0.95 


1.17 


1.19 




oo 


2.30 


4.96 


0.79 


0.14 


305 


0.54 


0.95 


1.15 


1.17 


5.309 




0.0 


1.56 


0.89 


0.04 


223 


1.0 


0.93 


0.98 


1.0 



at pn = p%- Since the equivalent scalar and vector potentials appearing in the wave equation 



for a point-like nucleon are Mjy — Mn and U v , respectively 



Eqs. (fy|) and 



imply 



large and canceling scalar and vector potentials for the nucleon in nuclear matter, 
potentials are comparable to those suggested by Dirac phenomenology 
calculations |L2| , and finite-density QCD sum rules 1J] , but smaller than those obtained in 



Such 
Brueckner 



QHD-I These potentials also imply a strong nucleon spin-orbit potential. Therefore, 
the essential features of relativistic nuclear phenomenology are recovered. The corresponding 
results for the nuclear matter compressibility, Ky 1 , are slightly larger than the corresponding 
value in the usual QMC model, but significantly smaller than that in QHD-I. The resulting 
total energy per nucleon for symmetric nuclear matter is shown in Fig. [I]. 

In the usual QMC model, the bag radius decreases slightly and the quark root-mean- 
square (RMS) radius increases slightly in nuclear matter with respect to their free-space 
values. When the bag constant drops relative to its free-space value, the bag pressure 
decreases and hence the bag radius increases in the medium. When the reduction of the bag 
constant is significant, the bag radius in saturated nuclear matter is 25 — 30% larger than 
its free-space value. The quark RMS radius also increases with density, with essentially the 
same rate as for the bag radius. This implies a "swollen" nucleon in nuclear medium, which 
has many attendant consequences ||T7| p3|| . It is also interesting to note that the result of 
25 — 30% increase in the nucleon size is comparable to those suggested in Refs. |]17|-|T9"| (see, 
however, 



Ref. pop. In the special case, g, 
in the nuclear matter is too large (B / Bq 



= and 6 = 4, we find that the reduction of B 
10%). This leads to unreasonably large values 



8 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 1. Energy per nucleoli for symmetric nuclear matter as a function of the medium density, 
with i?o = 0.6 fm and 5 = 8. Here the direct coupling model Eq. ( |3.lD for the in-medium bag 
constant is used. The solid curve corresponds to the usual QMC model, and the result from 
QHD-I is given by the open circles. The other three curves correspond to g%. = (long-dashed), 1 
(dot-dashed), and 2 (short-dashed), respectively. 

for the bag radius and the quark RMS radius in nuclear matter. 

The results corresponding to the limit of 5 — > oo [i.e., Eq. (|3.3|)1 are also listed in Table |. 
These results are not far from those with finite 5 values. As 5 increases, the results will 
approach saturation. The last row of Table [TI| gives the results obtained from fixing 
at its value predicted by the simple QMC model, 5.309 (for R = 0.6 fm). In this case, 
it is found that for any given 5 value the self-consistent solution requires = 0. When 
the coupling g% is tuned from zero to its corresponding value in the simple QMC model, 
the results interpolate between the QHD-I results and the usual QMC model results. If 
g% exceeds its value in the simple QMC model, the in-medium bag constant will increase 
instead of decrease relative to its free-space value, which is in contradiction with the physics 
discussed in the present paper. 

For a fixed g%, the coupling constants g„ and g u gets smaller as 5 gets larger. Recall 
that the reduction of B is controlled by both g^ and 5. While B / B and Mjy increase, U v , 
R/Ro, and Ky 1 decrease as 5 increases. From Table |, one can see that for a given value 
of B/Bq, one finds smaller Mjy and larger U v with a larger value of g q a . We also find that 
when 5 gets too small, a self-consistent solution no longer exists. 

For curiosity, we have also explored the results for negative 5 values. One can see from 
Eq. (|3.1|) that negative 5 values can also lead to a decreasing bag constant. In fact the 
physical quantities are continuous as 1/5 goes through zero. As 1/5 decreases from zero (5 
goes negative), the scalar and vector potentials continue the decrease seen in Table |. For 



9 



TABLE II. Coupling constants and nuclear matter results as obtained from the scaling model. 
The case of k = corresponds to the simple QMC model and the last row gives the result of 
QHD-I. Here the nuclear matter compressibility, Ky 1 , is given in unit of MeV, rj^ and r m denote 
the quark root-mean-square radius in nuclear matter and in free space, respectively, and U v = g^Zo 
is the vector mean field. The mass parameters are the same as in Table |. 



flo(fm) 
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^/4tt 




M* N /M N 


U Y /M N 




B/B 


x/xq 


R/Rq 


r ml r m 


0.6 





20.18 


1.56 


0.89 


0.04 


223 


1.0 


0.93 


0.98 


1.0 




1.0 


11.90 


2.27 


0.87 


0.06 


258 


0.87 


0.93 


1.02 


1.03 




2.0 


5.92 


3.60 


0.83 


0.10 


319 


0.69 


0.94 


1.08 


1.10 




2.95 


2.24 


7.78 


0.71 


0.22 


590 


0.36 


0.94 


1.27 


1.29 




3.0 


2.11 


8.32 


0.69 


0.24 


628 


0.33 


0.94 


1.30 


1.32 


0.8 





22.01 


1.14 


0.91 


0.03 


202 


1.0 


0.90 


0.99 


1.02 




1.0 


12.78 


1.76 


0.89 


0.05 


235 


0.89 


0.91 


1.02 


1.04 




2.0 


6.20 


2.88 


0.85 


0.08 


289 


0.73 


0.93 


1.08 


1.10 




3.0 


2.06 


6.39 


0.75 


0.18 


479 


0.42 


0.93 


1.24 


1.26 




3.1 


1.78 


7.32 


0.72 


0.21 


543 


0.36 


0.93 


1.28 


1.31 


1.0 





22.48 


0.96 


0.91 


0.03 


192 


1.0 


0.88 


1.0 


1.03 




1.0 


12.94 


1.54 


0.89 


0.04 


225 


0.89 


0.89 


1.03 


1.05 




2.0 


6.21 


2.58 


0.86 


0.07 


276 


0.75 


0.91 


1.07 


1.10 




3.0 


2.01 


5.68 


0.77 


0.16 


432 


0.46 


0.92 


1.21 


1.24 




3.17 


1.54 


7.18 


0.72 


0.20 


502 


0.36 


0.92 


1.29 


1.32 


QHD-I 




8.45 


12.84 


0.55 


0.37 


540 











g% = and 5 ~ —1.0, the resulting and U v are similar to those found in the usual QMC 
model. The compressibility is somewhat lower. 

All of the above results use the direct coupling model; we now present the results for the 
scaling model. The two coupling constants, g a and g u , are adjusted to fit the binding energy 
of saturated nuclear matter. There is only one free parameter, k, in this case. For various 
values of k and free space bag radius R , the resulting coupling constants and nuclear matter 
results are listed in Table |T[ 

Again, the decrease of the bag constant gives rise to the decrease of M^/Mjv and the 
increase of U v /Mn relative to their values in the simple QMC model. We see from Table 
1 1 that Mjy decreases and U v increases rapidly as k increases. For k = (2.95,3.10,3.17), 
corresponding to R = (0.6, 0.8, 1.0 fm), we find B/B ~ 0.36 at the saturation density. 
The corresponding results for and U v are 

M* N ~ 660 - 680 MeV , (4.3) 
U y ~ 190 - 225 MeV , (4.4) 

at pn = p%, which are similar to the direct coupling model results given in Eqs. Q4.1|) and 
(p~2"|), though the magnitudes for — and U v are slightly larger. These results are 
also consistent with those suggested by relativistic nuclear phenomenology and finite-density 
QCD sum rules. 
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FIG. 2. Energy per nucleoli for symmetric nuclear matter as a function of the medium density, 
with i?o = 0.6 fm. Here the scaling model Eq. ( |3.4| ) for the in-medium bag constant is adopted. 
The four curves correspond to k = (solid) , 1 (long-dashed) , 2 (dot-dashed) , and 3 (short-dashed) , 
respectively. The result from QHD-I is given by the open circles. 



One notices from Table [H] that the value of Ky increases quickly when k value is 
increased. This results from the increasing g w with increasing k. For the k values leading 
to Eqs. (|4.3| ) and Q4.4|) , the corresponding values for Ky are comparable to that obtained 
in QHD-I, which is too large compared with the empirical value. This feature is also shown 
in Fig. 0, where the total energy per nucleon for symmetric nuclear matter is plotted as a 
function of nuclear matter density for various k values, with R = 0.6 fm. The result from 
QHD-I is also plotted for comparison. The equation of state for the nuclear matter is much 
softer in the simple QMC model than in QHD-I. As k gets larger, the equation of state 
becomes stiffer. 

Both the bag radius and the quark RMS radius in nuclear matter are larger than their 
free-space values. As discussed above, this is due to the decrease of the bag constant in 



the nuclear medium. For the k values yielding Eqs. (O) and (|4.4|) , the corresponding bag 
radius and quark RMS radius at pn = p% are 25 — 30% larger than their free-space value. 
This is essentially the same as that found in the direct coupling model (see Table |). 

The sensitivity of our results to the free-space bag radius Ro is also illustrated in Table [T| 
For a given n value, the ratios B/B and M^/M N increase and the ratios R/Ro and U V /M N 
decrease as R increases. However, for the k values considered here, the sensitivity of our 
results to Ro is small. The sensitivity of the results from the direct coupling model to the 
choice of Ro is also small and similar to that in the scaling model. 
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V. DISCUSSION 



As stressed by Saito and Thomas H, in the simple QMC model, all the effects of the 
internal quark structure of the nucleon are summarized in the factor C(a). If C(a) = 1 is a 
constant, one would get exactly the same nuclear matter results as in QHD-I. In the simple 
QMC model, C(a) is much smaller than unity, which leads to much larger Mjy and much 
smaller U v than those required in QHD-I. 

In the present study, we introduce the medium modification for the bag constant. As 
shown in Eqs. (|3.2j ) and (|3.5| ), the effect of this modification is completely absorbed into the 
factor C(a). We observe that for various parameters considered here, both Eqs. ( |3.2j ) and 
( p.5|) lead to an increase in C(a). This indicates that the reduction of the bag constant in 
nuclear matter partially offsets the effect due to the internal quark structure of the nucleon. 
It is thus not surprising to find that our modified quark-meson coupling model gives smaller 
M^r and larger U v than those found in the simple QMC model. 

To further illustrate this point, we have plotted C(a) in Fig. [| as a function of g%a 
for the usual QMC model and for our modified QMC model with the scaling model for 
the in-medium bag constant. We see that in the simple QMC model, C(a) is small and 
decreases as g%a increases. The introduction of dropping bag constant gives an increase in 
C(a). When the reduction of bag constant is large, C(a) is approximately constant and 
significantly larger than one for small and moderate values of g^a; as gpr increases, C(a) 
decreases quickly. [The self-consistent solution requires g^a = (96, 82, 69, 58 MeV) at the 
saturation density, corresponding to k — (0, 1.0, 2.0, 3.0), respectively] A similar plot can 
also be done with the direct coupling model, but C(w) in this case is a function of g%, g^, 
and a, instead of g^a alone. 

We observe that the extent to which the bag constant drops in nuclear matter determines 
the physical outcome. Unless one expresses the bag constant in terms of QCD operators 
and solves QCD in nuclear matter, the change of the bag constant in nuclear medium is 
unknown. As such, one has to invoke model descriptions in order to obtain a quantitative 
estimate for the reduction of the bag constant in nuclear matter. 

In Ref. |]24|| , Adami and Brown have argued that the MIT bag constant is related to the 
energy associated with the chiral symmetry restoration (the vacuum energy difference be- 
tween the chiral-symmetry- restored vacuum inside and the broken phase outside). According 



to the scaling ansatz advocated by Brown and Rho |25| , the in-medium bag constant should 
scale like B/Bq ~ $ 4 , where $ denotes the universal scaling, $ ~ m*/m p ~ /*//„■•• 



which is density dependent. Here, the "starred" quantities refer to the corresponding in- 
medium quantities. This scaling behavior is argued to hold approximately at the mean-field 
level p^-p7|. Thus, one may get a rough estimate of the medium modification of the bag 



constant from the medium modifications of vector-meson masses which have been studied 
extensively |2l^pi2|. Taking the result for m*/m p from the most recent finite-density QCD 
sum-rule analysis ||29|| , we find $ ~ m*/m p ~ 0.78 at the saturation density, which gives 
rise to B/B ~ $ 4 ~ 0.36. This shows substantial reduction of the bag constant in nuclear 
matter relative to its free-space value. With this estimate, we obtain large and canceling 
scalar and vector potentials for the nucleon in nuclear matter, which are consistent with 
those suggested by relativistic nuclear phenomenology and finite-density QCD sum rules, 
though smaller than those found in QHD-I. This feature is seen in both models for the 
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FIG. 3. The factor C(w) as a function of g^W, with Rq = 0.6 fm. The solid curve is from 
the simple QMC model. The other three curves correspond to the scaling model, with k = 1 
(long-dashed), 2 (dot-dashed), and 3 (short-dashed), respectively. 



in-medium bag constant discussed here, implying a weak model dependence of our results. 

However, some caveats concerning the above estimate must be added. The Brown-Rho 
scaling is an ansatz based on the idea of partial chiral symmetry restoration in nuclear 
medium and the assumption that the scale anomaly of QCD could be modeled by a light 
dilaton field p5|]23 |. It is unclear whether this ansatz can be justified in QCD. Although it 
has been argued that many nuclear phenomena are connected to the partial chiral restoration 
in nuclear medium [2I|- |27 , 53| , 3"1 |, the only compelling evidence for partial chiral symmetry 
restoration in nuclear medium is that the magnitude of the chiral quark condensate, (qq) , is 

m 



substantially reduced relative to its vacuum value P5|-|55|. Recently, Birse has argued 
that the in- medium nucleon mass cannot be simply related to the change in the chiral 
quark condensate and there are other important contributions unrelated to partial chiral 
symmetry restoration. Moreover, it is known that both the MIT bag model and the QHD 
are not compatible with the chiral symmetry. 

Clearly, how the bag constant changes in nuclear matter is an important topic for further 
study. The investigation of finite nuclei and nuclear structure functions in the present 
model may offer some independent information and/or constraints on the modification of 
the bag constant. Work along this direction is in progress | 40fl . Another direction may be 
in motivating an explicit functional form for the in-medium bag constant, B(a), from a 
soliton-like model for the nucleon. Recently, an explicit form for B(a) has been suggested 
in Ref. |Tj|, based on the global color model of QCD. 



The quark-meson coupling model for nuclear matter is valid only if the nucleon bags do 
not overlap significantly. In his original paper H, Guichon has suggested Rq = 0.6 — 0.7 fm 
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FIG. 4. Quark wave functions as functions of radial coordinate. Here the scaling model Eq. (3.4) 
for the in-medium bag constant is adopted. The solid curves are the free-space quark wave 
functions. The other four set of curves correspond to, re = (long-dashed), 1 (dot-dashed), 2 
(short-dashed), and 3 (dotted), respectively. 



in order to keep the overlapping effect small. In our modified QMC model, both the bag 
radius and the quark RMS radius increase in medium relative to their values in free-space 
due to the dropping bag constant in medium. In Fig. f|, the spatial quark wave functions 
are plotted as functions of radial coordinate with the scaling model for the in-medium bag 
constant. We see that when the bag constant drops the quark wave functions are pushed 
outward. This depicts a "swollen" nucleon picture, which has important implications in 



many nuclear physics issues [17-23 



On the other hand, the increasing bag radius also implies larger overlapping effect than in 
the usual QMC model. When the bag constant in nuclear matter is significantly smaller than 
its free-space value (B/B ~ 0.4), we find R/Rq ~ 1.25 — 1.30 at the saturation density, 
which gives 4nR 3 p N /3 ~ 0.3 — 0.34 when Rq = 0.6 fm. This indicates that the overlap 
between the bags is still reasonably small at the saturation density, though a factor of two 
larger than in the usual QMC model. For larger Rq and/or higher densities, the overlap 
becomes more significant and the non-overlapping bag picture of the nuclear matter may 
become inadequate. However, it is unclear at this stage whether the overlap between the 
bags is already included effectively in the reduction of the bag constant and/or in the scalar 
and vector mean fields. Further study is needed to clarify such issue. We also note that in the 
scaling model with large k values, the resulting nuclear matter compressibility is too large 
compared to the empirical value. This may be fixed by introducing self-interactions of the 
scalar field, which, however, will introduce more free parameters. The direct coupling model, 
on the other hand, produces a more reasonable value for nuclear matter compressibility. 
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The QMC model is probably the simplest extension of QHD to incorporate explicit quark 
degrees of freedom, where the exchanging mesons are treated as classical fields in the mean- 
field approximation. To be more consistent, the explicit quark structures of the mesons 
should also be included and the physics beyond the mean-field approximation should be 
considered. It has been emphasized above that both the MIT bag model and the QHD are 
not chirally symmetric. As such, the quark-meson coupling models discussed here are not 
chiral models. At the hadron level, significant progress has been made in incorporating both 



chiral symmetry and broken scale invariance in relativistic hadronic models |41| , |42[| . So, the 
quark-meson coupling models may be extended by combining these hadronic chiral models 
and a chiral version of the bag model. Recently, it has also been argued that connections 
can be made between effective chiral Lagrangians and QHD model at the mean-field level 



VI. SUMMARY 

In this paper, we have modified the quark-meson coupling model by introducing medium 
modification of the bag constant. We proposed two models for the in-medium bag constant, 
direct coupling model and scaling model. The former couples the bag constant directly to 
the scalar mean field, and the latter uses a scaling ansatz which relates the in-medium bag 
constant to in-medium nucleon mass. Both models feature a decreasing bag constant with 
increasing density. 

The reduction of the bag constant in nuclear matter partially offsets the effect of the 
internal quark structure of the nucleon and effectively introduces a new source of attraction. 
This attraction needs to be compensated with additional vector field strength. The decrease 
of bag constant also implies the increase of bag radius in nuclear matter. This is consistent 
with the "swollen" nucleon picture discussed in the literature. 

When the bag constant is reduced significantly in nuclear matter with respect to its 
free-space value, we find that our modified quark-meson coupling model predicts large and 
canceling scalar and vector potentials for the nucleon in nuclear matter, which is qualitatively 
different from the prediction of the simple QMC model. These potentials are consistent with 
those suggested by the relativistic nuclear phenomenology and finite-density QCD sum rules. 
The internal quark structure of the nucleon seems to play only a relatively minor role. On 
the other hand, the reduction of the bag constant in nuclear medium relative to its free-space 
value may play an important role in low- and medium-energy nuclear physics. 
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APPENDIX: A 



In this appendix, we demonstrate that adopting Eq. with 5 = 4 for the in-medium 
bag constant and taking g%. = (and m° q = 0), one can reproduce the QHD-I results for 
nuclear matter. To this end, we show that the resulting expressions for the total energy per 
nucleon and the self-consistency condition for the scalar mean field are identical to those in 
QHD-I. 

In free space, the nucleon mass and the equilibrium condition for the bag are given by 
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bag 
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Solving these two equations, one finds that the combinations B Rq and M N R can be ex- 
pressed in terms of xq and Zq. This can also be seen easily from dimensional considerations. 
Since these two combinations are dimensionless, they must depend only on the dimensionless 
parameters xq and Zq. 

In the nuclear medium, the corresponding two equations become 
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Similarly, BR 4 and M^R are only dependent on x and Z. Since we take gl = and Z = Zq, 
x{= Xq) and Z are independent of the nuclear matter density. One therefore concludes that 
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Rq^ 4 
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Using Eq. (|3.1|) with 5 = 4, one gets 
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M N M N 

One can then rewrite the total energy per nucleon at finite nuclear matter density as 
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The self-consistency condition for the scalar field, Eq. (|2.13|) , can be expressed as 
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These two equations are identical to those required in QHD-I. Thus, fitting the nuclear 
matter binding energy at the saturation density, one should find the same scalar and vector 
couplings and hence the same strengths for scalar and vector mean fields as obtained in 
QHD-I. 
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